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Abstract. The subject of this article is a matched microstructure model 
for Newtonian fluid flows in fractured porous media. This is a homogenized 
model which takes the form of two coupled parabolic differential equations 
with boundary conditions in a given (two-scale) domain in Euclidean space. 
The main objective is to establish the local well-posedness in the strong 
sense of the flow. Two main settings are investigated: semi-linear systems 
with linear boundary conditions and semi-linear systems with nonlinear 
boundary conditions. With the help of analytic semigoups we establish 
local well-posedness and investigate the long-time behaviour of the solu- 
tions in the first case: we establish global existence and show that solutions 
converge to zero at an exponential rate. 



1 Introduction 

First ideas to use a two scale formulation to describe fluid flow in fractured porous 
media came up around 1960 e.g. by Barenblatt, Zheltov and Kochina ||9j. They reflect 
very well the exceptionell manner of the geometric conditions. The material possesses 
one natural porous structure and a second one is added by the dense system of cracks. 
A mathematical derivation of the model was later given in [5]. In the sequel weak 
formulations of the problem have been studied intensively either by homogenization 
theory [TSl [T^ or with the help of monotone operators [3 [TI] [2T] . Stationary solutions 
and the elliptic problem are for example treated in j50]. Reaction terms and evolving 
pore geometry have been considered in several papers [TH [THl [Ml US 123 [23 ■ Although 
the model in these studies has a similar form as in the present paper, the considered 
length scales are different. 

We consider the matched microstructure model (MM) as it was formulated by 
Showalter and Walkington in |29j . Assume we are given a macroscopic domain 17 C M", 
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and for each x G ft a cell domain ft^ C M". These ceh domains stand for the porous 
bfocks while fl contains in a homogenized sense the fissure system. 

The model (MM) consists of two parts: The macro model for a function u, that 
represent the density of the fluid on the domain f2: 

^^uit,x) - AMt,x) = f{t,u) + qiU)it,x), xen,teiO,T], 
u{t,x) = 0, xer,te{0,T], 

u(t = 0) = Uq. 



The micro model for the function U, that models the density in all blocks fla^ 



—u{t, X, z) ~ A.uit, x,z) = o, xen,z€n,,te{o,T], 



U{t,x,z) ^ u{x), X e fi,z eT^,t e {0,T], 

U{t = 0) = Uo. 

The coupling between the macro and the micro scale is reflected by two terms. Firstly 
the boimdary condition in the cells fix, 

U{x) = u{x) on dflx, for al\ x G ft (1) 

models the matching of the densities on the material interface. For this reason the 
model was introduced in |29] as matched microstructure model. Secondly the term 

qmt^x)^- ^EiL^ds^-^J^ Uit,x,z)dz. (2) 

represents the amount of fluid that is exchanged between the two structures. It acts 
as a source or sink term in the macroscopic system. 

Our interpretation of this model is based on the derivation of the coupled equations 
for the case of uniform cells at each point x in the considered domain ft. On this basis 
we will first present our restrictions on the geometry and the definition of suitable 
Banach spaces. Then we reformulate the problem as an abstract semilinear intial 
value problem on the product space Lp{ft) x Lp{ft, Lp{ftx)). Therefore we introduce an 
operator A that includes the highest order derivatives and the first coupling condition. 
In Theorem l3.4l we prove that — A generates an analytic semigroup which finally implies 
well-posedness of the matched microstructure problem. 

A further part of this work is the consideration of the long time behaviour of the 
solution. We show that for Dirichlct boundary condition the solution decays to zero 
at an exponential rate. In the Neumann case we prove mass conservation. Finally 
we consider a special two dimensional geometry and include a nonlinear boundary 
condition on dfl. A detailed derivation of this is given in [21]. We prove well posedness 
for weak solutions which can be improved using the concept of Banach scales. The 
approach is funded on work of Escher [13] and Amann [S] lU |6] . Our methods are quite 
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flexible and provide the high regularity of strong solutions. 

The outline of this paper is the following. First we present some basic lemmas 
for a uniform geometry. The largest section is then devoted to the variation of the 
cell domain and the transcription of the problem (MM) into an abstract initial value 
problem. There the mains result that —A generates an analytic semigroup is located. 
Further the spectrum of A is investigated for the case p = 2. The last chapter contains 
the model which includes nonlinear boundary conditions. With the help of a retraction 
from the boundary we transform the equations into an abstract semilinear problem, 
which leads to well-posedness of the original system. 



2 Some Aspects for Uniform Cells 

Let Q C M" be a bounded domain with smooth boundary F := and let (17.^, /i) 
be a measure space. Let X,Y be Banach spaces. If A is a closed linear operator 
from X to Y we denote with D{A) = {dom{A), || • the domain of definition of A 
equipped with the graph norm. Further we write U if we mean a representativ in Cp of 
a give function U G Lp. With [•] we indicate the equivalence class again. The shifted 
operators will always be denoted with bold letters. 

Lemma 2.1. Assume that {x A{x)) e C{Ti,jC{X,Y)). Let 

dom{A) = Lp{n,X), 

AU = \a{x)U{x) \ , forU G Lp{n, X), U €U. 

Then A is a well defined, bounded linear operator from Lp{D,,X) to Lp{n,Y). If 
further A{x) = A independent of x and A is a retraction, then A is a retraction as 
well. 

Proof. The continuity of {x t-^ A{x)) assures that A is well defined. Further we easily 
get AU e Lp{Q,, Y) and A is bounded because 

\\m\l^^n,x) = I \\A{x)U{x)ryd„{x) < maxP(x)||£ • \\U\\l^^^^^y 

Now assume! that A{x) = A is a retraction and let i? € L(Y,X) be a continuous right 
inverse of A, so that Ao R = idy. Let V € Lp{p.,Y), V S Lp(il, V) a representative 
of V . As before we define 

n €h{Lp{n,Y), Lp{n, X)), nv = [RV{x)]. 

Then a short calculation shows that this is a continuous right inverse for A. □ 

Lemma 2.2. For x € Q, let A(x) € A{X,Y) be a closed linear operator. Assume that 
there is Aq g A{X,Y), such that dom{A{x)) = dom{Ao), for all x £ fl. Furthermore 
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let {x ^ A{x)) € C(0,£(L»(Ao),y)). Then the operator 

dom{A) = Lp{il,dom{Ao)), 

AU = \a{x)U{x)\ , forU e dom{A), U e U. 

is a well defined, closed linear operator from Lp{D,, X) to Lp{fl, Y). If further dom{Ao) 
is dense in X , then A is densely defined. 

Proof. The first statement follows from Lemma 12.11 and a short calculation. The 
density assertion follows from the proof of Theorem 4.3.6 in [5]. □ 

Now we turn to sectorial operators. Let oj E R,9 E (0, tt). We set 

Se,u^ = {XeC;\^uj, \arg{\ - w)| < 9}. 

Lemma 2.3. Let the assumptions of the previous Lemma be fulfilled with X — Y . 
Assume further that there exist constants a; G M, G (0, tt), M > 1 such that for every 
X E fl, Se,ui C p{—A{x)) and 

||(A + A(x))-i||i(^) forXeSe,... 

Then A is sectorial in Lp{n,X). 

Proof Let A e Sq^^. With Lemma Owe define Ra G L{Lp{n,X)) by 

Ra{/:= [{\ + A{x))-^U{x)] 

for U E Lp{il,X). Then one easily calculates that this is the inverse of A + A. Thus 
A E p{—A). Furthermore it holds 

''~''~"il(x)I 

So we conclude that A is sectorial. □ 

3 The Semilinear Problem 
3.1 Geometry 

The main idea of this part is to relate the cell's shape to one standard cell, the unit 
ball B = B{0,1) C M". Let S = dB be its boundary. We assume that there are two 



||(A + A)-i[/||^^(^ ,,)< lj{X + Aix)r'\\l^^^\\U{x)r^df,{x) 

M 



n 
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mappings $ with 



$ : 17 X B ^ M" X M", 
{x,y) (a;, 

Now a cell at a point a; G is the image of B at x, i.e. il^; := ^{x, B). We set 

Q := IJ {a;} x VL^. 

Then Q = $(f2 x B). To assure that fix is a bounded smooth domain as well impose 
some properties of $, ^: 

^ eUp{nx B,Q), (3) 

<t>-^ eUp{Q,nx B), (4) 

•) e Diff(S,II:r), forahxef], (5) 

sup {\\d^'f{x)\\p,\\d:'f-\x)\\p} (6) 

2;Gf2,|a|<2 

Here || • ||p denotes the usual Lp-norm. The set Diff (S, il^;) shall denote all C°°- 
diffeomorphism from B to such that the restriction to the boundary S gives a 
diffeomorphism to Fa;. It follows from the assumptions that Q is measurable. Further 
for every x ^ fl, the set fl^ is a bounded domain with smooth boundary F^; :— dflx- 
Note that the special construction of $ implies that it is injectiv. Thus we will be able 
to work with the trace operator on B and transfer it to fix- The conditions ensure 
that the following maps are well defined isomorphisms. Given 2 < p < oo, we define 
pull back and push forward operators 

: Lp{n X B)^ Lp{Q) : U ^ U o 
$* : LpiQ) Lp{n X B) -.V o^. 

The following definition of a function space is based on Bochner's integration theory. 
In [T] (3.34, 3.35), it is proven for Sobolev-Slobodetski spaces that under these dif- 
feomorphisms Wp{B) is mapped onto Wp{ilx) for < s < 2. For s = we identify 
W;:(B) = Lp{B). The space Lp{yi^Wp{By) is now defined by means of the Bochner 
integration theory. We define 

Lp{n,w;{nx)) <^.{Lp{n,w;{B))). (?) 

We can prove that equipped with the induced norm 

ll/IU,. ll<I'*/l|L,(a,H^=(B)), / e Lp{^,W;{^x)). 
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this is a Banach spaceQ For the formulation of the boundary conditions on the cells 
we need suitable trace operators. Our assumptions ensure that we can restrict $*, <i>* 
to Lp(ri X S). We use the same notation for the puUback and push forward as on 
nx B. Let s > 0. We define 

Lp{n,w;ir,}) ■.= {Lp{n,w;{s))) 

\\U\\LAn,w^ir.)) = ||$*(7||i^(o,M/;(5)), U € Lp{n,W;iT,)). 

As before this is a Banach space. From Lemma 12.11 we deduce that the shifted trace 

tTs : Lpin, W^{B)) ^ Lpin, Wp'^iS)) : trs [/ = [trsU], 

is a well defined linear operator. The last trace in the brackets is the usual trace on 
B. Next we transport this operator to Q. We set 

tr : Lpin, W^in,)) ^ Lp{n, Wp'^r^)), tr := trs $*. 

The continuity of trs ensures that tr is a continuous operator. In particular 

trJ7 = implies trs{^*U) = 0. From Lemma \2A\ we conclude that iis is a retraction. 
There exists a continuous right inverse Rs of trs that maps constant functions on the 
boundary to constant functions on B. We define 

R := $*i?5$*. 

Then this is a continuous right inverse to tr. Let u E Lp{il). We identify 

Ru^R{u-ls)eLp{n,W^iB)). 

With Az we denote the Laplace operator in the coordinates z G ^l^- Similarly we write 
Ay and A^; for the Laplace acting on functions over _B or f2. The definitions above 
ensure that 

AzRu{x) = 0, for a.e. x(^n. (8) 

This will be helpful in later calculations. Another definition of function spaces for the 
matched microstructure problem can be found in [25) . 



Note that hypothesis l(6]l ensures that diflferent $'s within the class JSj to lead to equivalent 
norms. 
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3.2 Operators 

To use existing results for strongly elliptic operators, we first consider some auxiliary 
operators. Let Ai be the Dirichlet-Laplace operator on 57, 

dom{Ai) = w^{n) n w^^°{n), 

Aiu = — AjjM, for u G doni{Ai). 

It is well known that Ai is sectorial. For each x G J7 we define a Riemannian metric 
g{x) on the unit ball B. We write 

\/\g(x)\ := ydet5~(x), 

Then the regularity assumptions on $ imply that this metric is well defined and there 
exists constants Ci > 0,i = 1,2 such that Ci|a;p < J2i j 9^'' ^i^j — C'2|a;p. Let 
U € Lp{n, W^in^)), V = $*[/. We set 



qiU){x) := - / y/\g{^\g^={x)dy,Vix) ■ ds. (9) 



Here v = (z^i, . . . , Vn) denotes the outer normal vector on B. Then q(U) is a function 
in Lp(ri). Using the transformation rule for integrals one sees that this definition is 
consistent with ©. We define the operator A2 using the transformed setting 

dom{A2) = {ue Lp{n,w^{n,)y,tTU = o}, (lo) 

A2U = ^4A.,Vix)]. (11) 

The brackets [•] again indicate taking the equivalence class and V is a representative 
of V. Given x E fl, the operator Ax acts in the following way on w g Wp{B), 

•^^"^ = TfTrf (V\9{x)\g'Hx)dy^ V. 

V\9[x)\ ^ ' 

Note that Ax is the Laplace-Beltrami-operator with respect to the Riemannian metric 
g. It holds 

Lemma 3.1. The operator A2 is well defined. 

Proof. The coefficients of Ax depend continuously on x. Moreover the domain of 
definition is independent of x € 57. Since $ is defined up to the boundary of 51 the 
definition can be extended to its closure. So the hypothesis follows from Lemma [2.21 
and the properties of □ 

The following lemma collects some properties of the defined operators. Let R{X, A) — 
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{X + A) ^ denote the resolvent operator of —A for A G p{—A). 

Lemma 3.2. Assume that for any x ^ fl, '■= ^{x^-) is orientation preserving. 
Further assume that the Riemannian metric g^^ induced from $ is well defined. For 
each cell we define the transformation of the Dirichlet-Laplace operator , 

dom{Bx)^W^{B)f^Wl^°{B), 

B^v = AxV, for V e dom{Bx). 

It holds 

(a) The operators B^ are strongly elliptic in Lp(B). 

(h) Given x G fi, the operator B^ is sectorial. In addition there exists a sector 
Se,u, = {A e C; A 7^ \arg{\ - uj)\ < 9}, 
and a constant M2 > 0, both independent of x, such that 



p(-Bx) 3 S0^^, (12) 
M2 



l^(^>-S^)lli(Lp(i3)) < 7Y37-j: for all X e Sg^^. (13) 



(c) The operator B in Lp{fl x B), given by 

dom{B) = Lp{fl, W^{B) n Wl^\B)), 

BV=[BxV{x)], V e dom{B),V eV, 

is well defined and sectorial. 

(d) Set f := $*/ e Lp{Vl x B). If the function V £ Lp{n,W^{B) D W^^°{B)) is a 
solution of BV = /, then U := $»V^ fulfills 

-A,y{x, •) = fix, •), for a.e. x £ Q. 

Moreover 

U{x, z) — 0, for a.e. x G Q, z G T^. 

Proof. The first part foUows from the fact that strong ehipticity is preserved under 
transformation of coordinates. Part b) is a consequence of the definition of $ and 2.') . 
Theorem 3.1.3. With the help of Theorem 9.14, [T7] we conclude that the sector is 
independent of x. Now the rest follows by definition and Lemma [01 □ 

A more detailed proof can be found in [31. . Now we are ready to treat the coupled 
problem. Given u £ Lp{n), we set 

Doiu) := {U e Lp {n, W^in,)) ;tTU^u}. 
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This is a closed linear subspace of Lp{rt,Wp{ftx))- So we can define the operator A 
by 

dom(A) ~ [J {u} X Do{u), 

uew^{n)r]W^-°{n) 

A(m, U) = (^-A^u, [$,A**C/(a;)]) , for {u, U) G dom(A). 

Observe that the operator contains the matching condition ([T]). The exchange term 
q{U) will appear as a term on the right hand side of the abstract problem Let 
(/, 5) e ip(ri) X Lp(f2, Lp(17a;)), A e Sg^uj- We consider the system 

(A + A)(u,[/) = (/,g), for (w,C/) e dom(A). (14) 

This formally corresponds to 

XU~A,U^g, U€Dq{u). (16) 

Proposition 3.3. The operator ~A is the generator of an analytic semigroup on the 
space Lp(fl) x Lp{n, Lp{nx)). 

Proof. Let uji,0i such that Sg-^,i^-^ C p{~Ai), Se2.uj2 C p(— A2). Set 
= maxjoji, W2}, 6* = min{6'i, 02}- 

Then Sq^lo C p(— Ai) n/9(— A2). Take A G S'e,^;- Without restriction we suppose uj = Q. 
Since Ai is sectorial the function u = R{X,Ai)f solves ([T5|) . Furthermore there is 
Ml > 1 such that 

|A|||«||<|||A|i?(A,Ai)/||<Afi||/||. 

For U e i:'o(-«), it holds 

U - Rue Lp{n, Wpin^)) nkeiti = dom(A2). 

Here R is the extension operator defined in the previous chapter. So ^ implies that 
([TB]) is equivalent to 

\{U - Ru) + A2{U - Ru) = g - XRu. (17) 
Since B is sectorial, (jl7l) has the unique solution 

U = <I>,i?(A, B)$*(g - XRu) + 
So we have shown that ([H)) has a unique solution for A G Se^^;. Hence we conclude 
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that A G p(— A). To shorten the notation, we write Xq = Lp{il,, Lp{nx))- It holds 

\X\\\U\\ < |A|||$,i?(A,B)$*(7||xo + |A|||$,i?(A,B)$*Ai? R{X,Ai)f\\x„ 
+ \X\\\RR{X,A,)f\\x„ 
<Af2||**||||$1||l5llxo+M2||$,||||$*||||i?||Mi||/|U^(o)+Afi||i?||||/|^^ 

With this we estimate the norm of the resolvent of —A 

|A|||i?(A, A)|| = sup{|A|||C/|| + \X\\\u\\;u^ R{X,A,)f, 

U = R{X, A2)(5 - XRu) + Ru, \\f\\ + \\g\\ < 1}, 
< M2|l$*||||$*|| +MiM2||$*||||$*||||i?|| + +Mi =: M. 

Hence the sector is contained in the resolvent set p(— A) and the inequality above holds 
for some constant M > 1 independent of A. Hence A is sectorial and —A generates a 
holomorphic semigroup. □ 

We are now prepared to write the matched microstructure problem as an abstract 
evolution equation. Set w — (u, U), wo = (uo, Uq). We look for w satisfying 

(dtW + Aw^ f{w), te{0,T) ^^^^ 
\w{0) = Wo- 

To solve this semilinear problem we take ^ < 6 < 1. Our goal then is to show that 

/ : (0,T) X [Yo,D{A)]e ^ Yo := Lp{n) x Lp{n,Lp{n,)) 

is locally Holder continuous in t and locally Lipschitz continuous in w. For the initial 
value we will require that wq G [Yq, D{A)]q. Then results from Amann |P imply local 
existence and uniqueness. 

3.3 Interpolation and Existence Results 

Let Xo,Xi be two Banach spaces that form an interpolation couple. Let < O < 1. 
We denote with [Xq,Xi]q the complex interpolation space of order Q. Let trr : 

Wp{rt) — Wp " (r) be the trace operator on fl. R. Seeley showed in [2H] that 

He actually gives a proof for any normal boundary system (defined in the sense of [5S] , 
§3). To determine [Yq, D{A)]q we start with the case of uniform spherical cells B. Let 
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< e < 1. We set 

Xo=Lp{n)xLp{n,Lp{B)), 

Xi = {W'^iVL) n kertrr) x Lp{VL, W^iB) n kertrs). 

Then {Xq^Xi] is an interpolation couple. Due to Proposition 1.2.3.3 in [6 it suffices 
to interpolate both factors separately. Let 29 > ^. We deduce from (flQl) 

[Lp{^) ,W^{^)f^ ker trr] ^ = ker trr n W^^e 

Further the results in [10 and (flQl) show that 

[Lp{n,Lp{B)),Lp{n,W^{B)C^keYivs)\^^Lp[n, [Lp{B),W^(B)C^keYivs\Q). 

= Lp(f),%2Q(B))nkertrs. 

The map (/, <J>») : (u, U) i— > (u, $*(7) is an isomorphism. It maps 

Xo - Lp{n) X Lpin,LpiB)) ^ X =: Fq 

as well as 

Xi ^ (Wp2(^^) n kertrr) x Lp{n, W^{n^)) n kertr =: Yi. 
So {^0,^1} is an interpolation couple and Proposition I 2.3.2 from [B] implies 

for 20 > i. Finally we define the isomorphism 

J -.Yo^Yq: (u, U) ^{u,U + Ru). 

Here R is the retraction of the lifted trace. Then J maps Yi onto D{A). Clearly 
T : Yq ^ Yq : {u, U) 1— >■ {u,U — Ru) is the inverse of J. So J fulfills the conditions 
of Proposition I 2.3.2 from [B] for the interpolation couples {Fo;^i} and {Yq, D{A)}. 
Hence it maps [Yq, Fi]q onto [Yq, _D(A)]q . So for 26 > i, we get explicitely 

[yo,i?(A)]e = J([ro,yi]e)= (J M x {[/ e Lp(r!, W^^^e^^^)). ^ 

n kcr trr 

If 28 < ^ the boundary condition drops in both scales. Hence we conclude 

[YQ,DiA)]^ = W^'^iQ) X Lp{n,w^''{n,)). 
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A similar analysis can be done for real interpolation functors. In particular 

(Fo,i?(A))i_i_p - U {u}x{U eLp{^,wl'~^{^^)yMU ^u}. (20) 

ueWp " (O) 
n kcr trr 

Let < 9 < 1 and write X'^ ~ [Yq, D{A)]q. We consider functions 

/ = (/,.g):[0,oo)xX®->yo- (21) 
Theorem 3.4. Let f — {f,g) be as in (I?!]) . Assume 

/eci-([o,^)xx«,yo), 

is locally Lipschitz continuous for some < < 1. Then for any (uq, Uq) G , there 
exists T = T(uo,C/o:0) > 0; such that \18\) has a unique strong solution w = {u,U) 
on (0,T) which satisfies initial conditions 

u{t = 0) = uo and U{t = 0) Uq. 

In particular 

u;eCM([0,T),Fo)nC([0,T),X®). 

Proof. With the above considerations, Theorem 12.1 and Remark 12.2. (b) from [S] 
can be applied to the abstract equation (1181) . The regularity results are proved in 

m- □ 



In order to state our main result on the matched microstructurc model let 

Xf ~ %'®(r!) n kertrao, ® > ^' 

denote the first component of the interpolation space X^. 
Corollary 3.5. Let 



fee'- ([o,oo)xxy^+^/^i,(i])) 



be given. Then for (uq,Uo) G X^/^+^/^p, there exists T > 0, such that the matched 
microstructurc problem (MM) has a unique .strong solution on (0,T). 

Proof Let e = i + ^. Let U G Lp{Q,W^^{n^)), V = Since q{-) is time 

independent it remains to show the Lipschitz continuity in U. There is C > 0, such 
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that 

< 



p 

dx 



yWx^WHx, s)dy^V{x, s)vj ds 
<cl ma_x _ {^W^\g^={x,y)\Y [ J2 f \dyM^,s)\'dsdx 

i, j K — ± 



< C 

The constant Cp is the embeddmg constant of Lp{B) into Li{B). Together with the 
hnearity of q this shows that q is locahy Lipschitz continuous in 

Wp^^n) X Lp{n,Wp^^n^)) d x^. 

Setting finahy f{u,U) := {f{u) + q{U),0), the assumption fohows from Theorem 

[231 ~ □ 

3.4 Exponential Decay under Dirichlet Boundary Conditions, Neumann BC 

Now we assume that there are no external sources in the system, i.e. / = 0. We 
will show that the corresponding solutions decay exponentially fast to zero. First we 
investigate the spectrum of A. Let p — 2. The space L2{^) x L2{^, L2{B)) is a Hilbert 
space. This implies that Yq is a Hilbert space with the inner product 

((u, U), {w, W))y, = (u, w)L,in) + i'^*U, $,W^)L,(j7xi3) 

for (u, U), {w, W) e Yq. We introduce an extended operator by 

dom{Aq) = (iom(A), 

A,(m, U) = A{u, U) - {q{U), 0), for aU {u, U) G dom(AJ. 

We investigate the spectrum of Ag. It is convenient to introduce a weighted space Yg. 
Set 

Yg = X L2{n, L2{n,, VWD), 
llKt/)ll?', = lkllL(o) + ll'i>*C/|lL,,,,(,,^,. 

This is a well defined Hilbert space with respect to the inner product 



{{u,U),{w,W))y, = UW+ y^\g\^*U<i>*V, for {u,U), {w,W) e Y^ 

Jn JnxB 

We can show that with this modified inner product Yg is also a Hilbert space. 
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Lemma 3.6. The operator Aq is self adjoint in Yg. 

Proof. Take (u, U), {w, W) G dom(A,). Let V = ^*U, Z = ^*W. Then 

and 

{Ag{u,U),{w,W))Y,^ I {~^Mx)w{x)~q{U){x)w{x)) dx 
Jn 




Manipulation of the last integral in (12^ by partial integration together with the bound- 
ary conditions on fl and B shows 




= - / / V{x,y)J2dyJg''V\9\dy.Z(x,y)) dydx+ f q{U)w ~ [ u q{W). 
Jn Jb , w J^ Jit 

This implies that is symmetric. From the theory of elliptic operators and the 
representation of A we can conclude that A is invertible, in particular we have that 
imA — Y. Our next step is to show that this holds for A^ as well. Let {v, V) E Y. 
We show that there exist {z,Z) G dom{Aq) with Ag{z,Z) = {v,V). First we know 
that there are {u,U) G dom{A) such that A{u,U) = {v,V). Then 

Aqiu,U)^{v-qiU),V). 

Clearly {q{U),0) G Y. Again there are functions {w,W) £ dom{A) with A(w, W) = 
(g(C/),0). This implies that W{x) = const, for a.e. x E VI. Thus q{W) = 0. Then 
from the linearity of Aq it follows that 

Aq{{w, W) + {u, U)) = {q(U), 0) + [v~ q{U), V) = {v, V). 

Thus Aq is symmetric and im{Aq) = Y . It follows from the fact that 

ker(A;) = im{Aq)^ = {0}, 

that the dual operator is injektiv. Thus Aq C A* implies the assertion. □ 

Lemma 3.7. It exists a constant a > such that 

i-Aq{u, U), {u, U))y, < ~a{{u, [/), {u, U))y^ , for all {u, U) G dom{Aq). 

Proof. Let {u, U) G dom{Aq). Set V = $*[/. We make use of equivalent norms in 
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Sobolev embedding results and that the metric g^^ is bounded from below. Let C > 
denote an appropriate constant. It holds 

{~A,{u,U),{u,U))y^ 

Jn Jn Jn Jb ■ ■ 

''■J 

—u{x) 

J29''V\9\dy.V-i^,V{xj 



= ~q{U) 

< -C + niL(aL.(a.))) = -'^ (K U), {u, U)) . 

Hence we have obtained a bound for the numerical range of — Ag in the weighted 
space. □ 

The spectrum of a self adjoint operator is contained in the closure of its numerical 
range (see [35], Section V, §3). Hence the spectrum of — lies totally on the right 
hand side of — tr. Since the weighted norm and the usual norm on Y are equivalent, 
we also get a spectral bound for — Ag in the unweighted space. So the right half space 
is containt in the resolvent of — Ag. We set 

Q : D{Ai) ^ Yo : {u, U) ^ {-q{U), 0). 
Then Q e C{Yi,Yo). Obviously 

-Ag = -A + Q 

and the conditions of Proposition 2.4.1 in [33] are satisfied. So Ag is sectorial. The 
matched microstructure problem is equivalent to 

idt(u,U)+A,{u,U)^0, te(0,T), 
\iu,Um^{uo,Uo). 

Proposition 3.8. Let {u, U) he a solution of ^MW- 
Then (u, U) \, (0, 0) exponentially fast. 

Proof. This follows from the the fact that for analytic semigroups the growth bound 
and the spectral bound coincide as it is e.g. shown in |12^ Corollary 3.12. □ 



Proposition 13.81 allows also to apply the principle of linearized stability to the semi- 
linear version of (MM), provided / is of class C^, cf. [33]. We can also treat a 
modified model with no-flux or Neumann boundary conditions on fl. Thus we want 
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d^u = on r. We set 

dom(Af ) = {u e W^{VL);d^,u = on r} , 

A^u = -I^xU, for all u E dom{A^). 

The boundary conditions in the cells are not changed. Hornung and Jager also derived 
this model in i20j . Then —A^ is the generator of a strongly continuous, analytic 
semigroup in Lp{rt) for 1 < p < oo. We set 

dom(A^) = U ^ ^o(")' 

A^{u, U) = {A^u, [^^AxV{x)]), for (u, U) G dom(A^), F = 

The modified model can be formulated as evolution equation 

dt{u,U) + A''iu,U) = {q{U),0), iG(0,r), 
(u,t/)(0) = (uo,C/o). 

The changes in the operator occur only on the macroscales. Thus they can be treated 
with well known results for elliptic operators on bounded domains. So the same con- 
siderations as for A can be done for A^. Existence and uniqueness can be proved 
similarly as in Chapter 13.31 Nevertheless the qualitative behaviour is different. 

Proposition 3.9. Let (uq, Uq) £ Wp{n)xLp{Q, Wp{Qx)) and let {u, U) be the solution 
to the matched micro structure problem with Neumann boundary conditions {24^ on 
some time interval [0,T]. Then the material value 

S{u,U):= f u+ f f ^\^*u, ie(o,r], 

Jn JnJB 

is preserved. 

Proof. This follows from a straight forward calculation □ 



4 Nonlinear Boundary Conditions 

The following version of the MMP brings gravity into the scheme. Let us restrict 
ourselves to uniform cells. Assume we are given a function / : (0,27r) (0,oo) 
periodic and differentiable. We consider the fixed domain 

nf = {ix,y)e§'xR-0<y<f{x)}. 

It is shown in Figure [T] The gravitational force points into the —y direction. The 
almost cylindrical domain Vlf can be treated with the same methods as before. For a 
work on the torus see [14]. Let /i : 17/ x (0, T) — M describe the sources and sinks in 
the macro system. A solution of the matched microstructure model with gravity is a 
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f(x) 



Figure 1: The periodic domain flf 



pair of functions (u, U) that satisfies 
dtu 



(P) 



dtU 



u 

AyU 

U 

iu,um 



h + q{U), on G (0,r), 

-u^ onro,t e (o,r), 

Po, onr/,i e (o,r), 

0, in flf X B,t e {0,T), 

u, on nj X S,t e {0,T), 

{uq,Uo), on rif X (flf X B). 



Let V = u - po ■ trtf , V = U - pot 
q{V + po) = q{V). So {v, V) solves 



nfxB- 



By definition © it holds that q{U) = 



(P') 



dfV — AxV 

V 

dtV - AyV 

V 



h + q{V), 
-{v + p^f, 
0, 
0, 

("0 - Pa, Uq - Pa), 



on nf,te (0,T), 
on ro,te (0,T), 
onr/,i e (0,T), 
in fi/ X e (0,T), 
on rif X S,te (0,T), 
on Qf X {flf X B). 



To treat the nonlinear boundary condition in the macroscopic scale we us a weak 
formulation. As before we define operators Ai and Ag with linear zero boundary 
conditions {d2V = on Fq, v = on F/). The mixed Dirichlet-Neumann conditions 
on ri/ do not effect the properties of the operators. Especially Aq is selfadjoint and 
— Ag is the generator of an analytic semigroup in the Hilbert space 

Yo=L2{nf) X L2{nf X B). 

The operator Ai is a well known form of the Laplace operator and so it is invertible. 
The presented method is due to Amann [3] and Escher [T3] . The main idea is to move 
the nonlinearity from the boundary to the right hand side h. Therefore we need to 
construct an appropriate inverse operator from the trace space on Fq to the domain Q, f . 
The resulting semilinear evolution equation can be treated as before. For u € L2{^f) 
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we set 

D'o{u) = {Ue L2inf,Wi'iB)),tTs U = u}, i < s < oo. 

In the rest of the chapter we drop the index / of 51/. Let tro,tr/ denote the trace 
operators onto Fq and Tf. With the boundary operator on we mean an operator B, 
with Bu = tro di,u + tifu. We define 

( {iu,U),ue Wi%n),U (z D^{u),Bu^O} , for I < 2s < oo, 
1;= < {(u,t/),ue tr/M = 0}, fori<2s<|, (25) 

[H/|"(r2) X L2{n,Wi'{B)), for < 2s < i. 

To construct a suitable retract, we first restrict ourselves to the macro scale. Let 
Ai ~ —Ax- Considered as an unbounded operator in L2{i^) it is closable. Together 
with B it fits into the scheme of [13], Chapter 3. We will use the same notation. Let Ai 

be the closure of Ai. Then VF|(ri) A- D{Ai). In addition we set Cu = tro it + tr/ d^u, 
and 

dWi' = 1^2^'"^ (Fo) X W^'^hTf), diWl' = T4^f "^(Fo) x wl'^^iVf), 

for < s < 1. Combined together, the map {B,C) G lj{W^{n),dWi x diW^) is a 
retraction. So we can apply Theorem 4.1. from [3]: 

Proposition 4.1. There exists a unique extension ('B, C) G L{D {Ai) , SVF" ^ SiIVj) 
of {B,C) such that for u G D{Ai), v G VF2^(r2) the generalized Green's formula 

{v,Aiu)yo + {Cv,Bu)qwO = {AiV,u)yo + {Bv,Cu)q^wo 

is valid. 

Then from interpolation theory (Proposition I 2.3.2 in [6 ) and well known a priori 
estimates for Ai, it follows that 

(Ai,B) G Isom{{D(Ai),Wi{n))g,L2{n) x 9 G [0,1]. 

Therefore we can define the right inverse 

Re - (Ai^By^ |{0} X aw^l^ 9 G [0, 1]. 

Then Rg G L{dW2^ , VF|*(r2)). We now add the microscopic scale. With 3^" we mean 



{{u, U) G Wi'{n) X L2{n, Wi'iB)); U G i^g^(u)}, i < 2s < 2, 
Wi'{VL) X L2(ll, W^^iB)), < 2s < i. 



So if s < I the two sets Ys and 3^" coincide. Define Hg G h{dWi^ , y'^^ ) by 

Kgu^ {RgU,RgU-lB), U & dWf . 
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We set 

doW^'' ^{ue dWi'^; tTfu = 0}. 
Obviously it is a closed linear subspace of 9PV^|^. It can be identified with the space 
W^^'^Tq). So it holds 

Iie{doWf )cY2e, if 20 < | 

For the formulation of the abstract evolution problem we use the scale of interpolation 
and extrapolation spaces {(Yq, {Aq)a),a G M} as it was defined by Amann [3]. For 
< a < 1 this corresponds to the interpolation spaces in Section 13.31 and Definition 
([251) . We set as in [13] 

A-(A,)„i, H = Y_., D = Y.=DiA). 
Then the duality theory tells us that D ~ H' and the duality pairings satisfy 

{u,if}H = {u, v)yo = {u\v)yoi for u e D,v (EYq. 
Let a : D X D ^ M. he the coercive bilinear form 

a{u,v) = / \/xU-\/xvdx+ / \/ zU ■ \7 d{x, z), il,v d D. 
Jn JnxB 

Usually it is clear from the context whether we refer to the function u living on or 
the pair u = (u, U). For u,v d D we get 

{v, Au)h = I V.,v- V^u - f vq{U) + f V,V- - f VV,U ■ v 
Jn Jn JnxB JnxS 

— / vdijU— / v\/ xU ■ V — a[uTv). 

Jvo Jvf 

The possible approximation of u by functions in Yi and the continuity of the left 
and right hand side justify this formal calculation. To treat the nonlinear boundary 
condition we define the map G : Z? — > L2{Y'q), 

G{u) = -tro(^t + po)^• 
We have to show that (/i, G) satisfies the assumption (3.6) of [13]. For h we just assume 
that /i e lo- For G the properties are summarized in the following lemma. 

Lemma 4.2. 

for any fixed /3 G (— ^, — ^),, and the Lipschitz continuity is uniform on bounded sets. 

Proof Fix /3 e (-5,-3)- Let u e D. Then also u + po e W^{n). From g], 
Theorem 4.1 and the fact that the Besov space i?|2(^^) = 1^2(^)1 know that the 
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multiplication W2^(r2) • W2i^) — J> W2 ^(^1) is continuous for < e < 1. We conclude 
that for fixed e < ^ we have 

{u + pof e Wt'i^) and tro(M + pof G wj~'{To). 

Then by Sobolev embedding it holds wj^^'iVo) A L2(ro) A T4^|'^+^ (Fo). The second 
inclusion follows from the definition of W2 as a dual spaces for s > 0. So finally 

-tioiu + po'f e w^f +^(ro). 

The Frechet derivative of G is the linear operator dG{u)v = — 2tro(w + po)v. Thus 

G e Ci(I?, VK2'''^^(ro)). It remains to show that the map is uniformly Lipschitz 
continuous on bounded sets. Let W G D he bounded. Take u,v ^W. Then 

tro(u + po)^ - tro(u + pof = tro - tro + 2po(tro u - tro v). 

Clearly the last term is uniformly Lipschitz on W . Further W2^(ri) ^ C{Vl). So a 
bounded set in D is bounded in C(ri). Thus there exists a constant ci > 0, such that 
||w||oc < ci for all u e M^. It follows from this and Sobolev embeddings that 

II tro M^- tro ?;^|| 2/3+1 ^ < C|| tro(M^ - w^)||L,(ro) 

(lo) 

< L\\u-v\\^ri(n)- 

The last constant L is independent of m, w G W^. This completes the proof. □ 

Now we define the right hand side to write (P') as an abstract evolution equation. 
Let R := Ri . Then it holds for u € 1^2^ (Lq) and v <E D 

{v,ARu)h = {v,u) 1 =: {v,u)ro, (26) 
W2 (ro) 

in the sense of trace. We set 

= (/i,0) + ARG(u), ueD. 

Note that the second component of F{u) vanishes since R^u ■ Is is constant on each 
cell. By assumption h G Yq ^ Yp. It was shown in [2], p.301, that under this 
circumstances F is well defined and the previous lemma ensures that 

F G C^{D,Yp) is uniformly Lipschitz continuous on bounded sets. 

Proposition 4.3. For each uq d D there is a unique maximal solution u(-,uo) G 
C{[0,Ti),D) of the semilinear Cauchy problem 

'il + Au = F{u), u(0) = ito (27) 
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with < Ti < CX3. In addition for any e € (0, i) holds 

u e c((o, Ti), 1 ) n ci((o, Ti), 1 ). 

Proof. Take /? = —^+e. Then the assertion follow from [3], Sect. 12 and the previous 
lemma. □ 

By a weak solution of (P') we mean a function u G C^([0, Ti), Z?) such that the 
initialcondition u{0) = uq — p is satisfied, and 

{iP,u)h + a{ip,u)dt ^ / ( ((^, (/i, 0))h + / (pG{u)] dt + {Lp{0),uo - po) h 



for aU < T < Ti, e C([0,r],i:>) n C^{[Q,T],H) with (^(T) = 0. So the above 
considerations show that 

Corollary 4.4. 

For each Hq G Yj. i/iere exists a unique maximal weak solution of (P')- 
Proof. This follows from the representation of A and (PSI) . □ 
Remarks: 

(i) The construction in [TB] allows to consider a more general h E C^(-D,y^). This 
means a full semilinear version of (P') can be treated. 

(ii) Abstract results on evolution equations in interpolation-extrapolation scales ensure 
that the solutions satisfies 

(u, U) e C((0, Ti), Fi) n ^1(0, Ti), Fo). 

So the system (P') is satisfied pointwise in time. 
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